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^H ' Abstract: 



The process pp — > ttbb + X represents a very important background reaction to 
searches at the LHC, in particular to ttH production where the Higgs boson decays into 
a bb pair. A successful analysis of ttH at the LHC requires the knowledge of direct ttbb 
production at next-to-leading order in QCD. We take the first step in this direction upon 
calculating the next-to-leading-order QCD corrections to the subprocess initiated by qq 
annihilation. We devote an appendix to the general issue of rational terms resulting from 
ultraviolet or infrared (soft or coUinear) singularities within dimensional regularization. 
There we show that, for arbitrary processes, in the Feynman gauge, rational terms of 
infrared origin cancel in truncated one-loop diagrams and result only from trivial self- 
energy corrections. 
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1 Introduction 

The search for new particles will be the primary task of the LHC experiment at CERN 
starting this year. The discovery of new particles in the first place requires to establish 
excess of events over background. The situation at the LHC is particularly complicated 
by the fact that for many expected signals the corresponding background cannot entirely 
be determined from data, but has to be assessed upon combining measurements in signal- 
free regions with theory-driven extrapolations. To this end, a precise prediction for the 
background is necessary, in particular including next-to-leading-order (NLO) corrections 
in QCD. Since many of these background processes involve three, four, or even more 
particles in the final state, this kind of background control requires NLO calculations 
at the technical frontier. This problem lead to the creation of an "experimenters' wish- 
list for NLO calculations" at the Les Houches workshop 2005 [1], updated in 2007 [2], 
which triggered great theoretical progress in recent years (see for instance Refs. [1-15] 
and references therein). Meanwhile the listed processes involving at most five particles in 
loops have been completed in NLO QCD, including the production of WW -|- jet [16, 17], 
weak-boson pairs plus two jets via vector-boson fusion [18], and triple weak-boson pro- 
duction [9, 19]. However, none of the true 2 — >• 4 processes has yet been addressed at 
NLOlj. Among those processes, pp -^ ttbb + X has top priority. This process has also 
been discussed as signal of strong electroweak symmetry breaking [21]. 

The process of ttbb production represents a very important background to ttH produc- 
tion where the Higgs boson decays into a bb pair. While early studies of ttH production 
at ATLAS [22] and CMS [23] suggested even discovery potential of this process for a light 
Higgs boson, more recent analyses [24, 25] with more realistic background assessments 
show that the signal significance is jeopardized if the background from ttbb and tt + jets 
final states is not controlled very well. This is a clear call for improved signal and back- 
ground studies based on NLO predictions to these complicated processes. For the ttH 
signal [26,27] and the tt + Ijet background [28] at the LHC such predictions have been 
accomplished in recent years. 

The dominant mechanism to produce ttbb final states in hadronic collisions is pure 
QCD. In leading order (LO) quark-antiquark (gg) and gluon-gluon (gg) initial states 
contribute, where the latter strongly dominate at the LHC because of the high gluon 
flux. Being of order a^ the corresponding cross sections are affected by a very large scale 
uncertainty, which amounts to a factor two or more. Technically the qq channel is simpler 
to deal with — though still demanding — and thus represents a natural first step towards a 
full treatment of pp -^ ttbb + X at NLO. In this paper we report on this first step and 
present some details of the calculation as well as numerical results. These results do not 
yet significantly improve the predictions for the LHC, but on the one hand form a building 
block of the full calculation and can serve as benchmark results for other groups on the 
other. Moreover, this step proves the performance of the applied strategy and methods, 
providing confidence that the more complicated gg channel can be attacked widely in the 
same way. 

In Section [2] we give a brief description of the NLO calculation, followed by numerical 
results on integrated cross sections in Section [31 Appendix [X] provides a general discussion 



^Progress in the calculation of the virtual corrections to uu -^ ssbb was reported in Ref. [20]. 



of rational terms in one- loop amplitudes, and Appendix [B] outlines some technical details 
concerning our treatment of the Dirac algebra. Finally, as a benchmark, in Appendix O 
we give numerical results for the matrix element squared in lowest order and including 
virtual corrections for one phase-space point. 

2 Details of the calculation 

In LO QCD seven different Feynman diagrams contribute to the partonic process 
qq — > ttbb; the various topologies are shown in Figure [TJ The virtual QCD corrections 
comprise about 200 one-loop diagrams, the most complicated being the 8 hexagons and 
24 pentagons, which are illustrated in Figure [2J The real QCD corrections in the qq 
channel are induced by gluon bremsstrahlung, qq —>■ ttbbg, where the corresponding 
64 diagrams are obtained from the LO graphs upon adding an external gluon in all possible 
ways. In the following we briefly describe the calculation of the virtual and real corrections, 
where each of these contributions has been worked out twice and independently, resulting 
in two completely independent computer codes. 

2.1 Virtual corrections 

The general strategy for the evaluation of the one-loop corrections is based on the 
reduction of the amplitude A^*^^^ of each (sub)diagram F in the following way. 
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where the colour structure C*^^) present in the (sub)diagram is factorized from the remain- 
ing colour-independent part. The decomposition of the colour structure. 
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is done in a basis {Ck} consisting of six elements, which can be chosen as 
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Figure 1: Two different diagram topologies contributing to qq -^ ttbb in LO QCD; there 
are six explicit diagrams of the first and one of the second kind. 





Figure 2: Diagram topologies for pentagon and hexagon graphs contributing to qq 
at one loop in QCD; there are 24 explicit pentagons and 8 hexagons. 
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Here T'*, f"'^'^, d°''"^ are the usual SU(3) objects, and the tensor products connect the three 
fermionic chains. The LO amplitude is decomposed in colour space as 
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In general each loop diagram gives rise to 3"" colour- factorized amplitudes of type (12.11) . 
where n4 is the number of quartic gluon vertices in the diagram. However, for most 
diagrams n4 = 0, and the colour structure factorizes completely. The colour separation 
implies that the computation time for individual loop diagrams does not scale with the 
number of colour structures present in the basis {Ck}. 

The colour- free parts of A^^'"^ are written as a linear combination of so-called standard 
matrix elements (SMEs) J^m,{{Pa}), which contain all Dirac chains and the polarization 
information. Since the computing time scales with the number of SMEs, it is important to 
reduce the set of SMEs {Aim} as much as possible. To this end, we employ an algebraic 
procedure based on four-dimensional relations that are derived from Chisholm's identity 
whenever their use is admitted. For massless external fermions this four- dimensional 
reduction has been described in detail in Sections 3.1 and 3.3 of Ref. [29]; here we had 
to generalize this approach to one massive and two massless spinor chains. In this way 
some thousand different spinor chains are reduced to about 150 SMEs M.m{{Pa})- A brief 
description of this procedure, which is implemented in two independent Mathematica 
programs, is outlined in App. [Bl 

The one-loop correction to the spin- and colour-summed squared amplitude induced 
by r reads 
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where the interference of the LO amplitude with the elements of the SME and colour 
basis. 
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has to be calculated only once per phase-space point. Moreover, the colour-correlation 
matrix X^coi ^k^i obviously is independent of the kinematics and is only calculated once 
and for all. The most time-consuming components of the numerical calculation are the 
scalar form factors J-'^\ which are linear combinations of the Lorentz- invariant coefficients 
of A^-point tensor loop integrals with rank R < 3 and degree A^ < 6, 

^in\{{PaPb)}) = j: E J^Sn,...,jMP^P'^)})TL.,MP'^P'^)})- (2-7) 
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The evaluation of one-loop tensor integrals Tj^ ■ follows the strategy of Refs. [4, 5p that 
was already successfully used to compute the NLO electroweak corrections to e"'"e" -^ 
4fermions [29,31]. In this approach the analytic expressions are not reduced to master 
integrals. In contrast, the tensor integrals are evaluated by means of algorithms that 
perform a recursive reduction to master integrals in numerical form. This avoids huge 
analytic expressions and permits to adapt the reduction strategy to the specific numerical 
problems that appear in different phase-space regions. 

The scalar master integrals are evaluated using the methods and results of 
Refs. [32,33]. Ultraviolet (UV) divergences are regularized dimensionally in both eval- 
uations, but infrared (IR) divergences are treated in different ways as described below. 
Following ideas from the 1960's [34], tensor and scalar 6-/5-point functions are directly 
expressed in terms of 5-/4-point integrals [4, 5]|f| Tensor 4-point and 3-point integrals 
are reduced to scalar integrals with the Passarino-Veltman algorithm [35] as long as no 
small Gram determinant appears in the reduction. If small Gram determinants occur, two 
alternative schemes are applied [5]jj One method makes use of expansions of the tensor co- 
efficients about the limit of vanishing Gram determinants and possibly other kinematical 
determinants. In the second (alternative) method we evaluate a specific tensor coeffi- 
cient, the integrand of which is logarithmic in Feynman parametrization, by numerical 
integration. Then the remaining coefficients as well as the standard scalar integral are 
algebraically derived from this coefficient. The results of the two different codes, based on 
the different methods described above are in good numerical agreement. Although both 
versions of the virtual corrections basically follow the same strategy for the evaluation 
of loop integrals, they are based on independent in-house libraries. In each of the two 
calculations the cancellation of IR and UV singularities was checked with high precision 
in the numerical results. 

Version 1 of the virtual corrections starts with the generation of Feynman diagrams us- 
ing FeynArts 1.0 [36]. Their algebraic reduction is completely performed with in-house 
Mathematica routines. In detail, D-dimensional identities (Dirac algebra, Dirac equa- 
tion) are used until UV divergences cancel against counterterms. IR (soft and coUinear) 
divergences are regularized dimensionally and separated from full diagrams in terms of 
3-point subdiagrams as described in Ref. [37]. We subtract the IR-divergent part A^^ing 
from the amplitude Ad^^'^'^ of a (sub-)diagram F, where D indicates dimensional regular- 



^We note in passing that the reduction methods of Refs. [4,5] have also been used in the related 
calculation [30] of NLO QCD corrections to the 2^4 particle process 77 — > ttbb at a 77 collider. 

■^Similar reductions are described in Ref. [7]. 

"^Similar procedures based on numerical evaluations of specific one-loop integrals [3, 7] or expansions 
in small determinants [6] have also been proposed by other authors. 



ization, and add it back. Note that M-shlg can be easily constructed already at the inte- 
grand level of the whole diagram following Ref. [37]. In the IR- finite and regularization- 
scheme-independent difference Ai^^'^^ — TVfgjng we can then switch from dimensional 
regularization to a four- dimensional scheme, 

^(r,^) = (^(r,D) _ _^(r.D)) ^ _^(r,D) ^ (^(r,.) _ _^(r,A)) ^ _^(r,z.)^ (2.8) 

where A indicates any mass regulators in four dimensions. Specifically, we introduce 
infinitesimal light-quark and gluon masses with the hierarchy nig <^ rriq to regularize the 
IR singularities. The evaluation of Ai^^'^^ then proceeds in D = A — 2e dimension merely 
to regularize the UV singularities, i.e. so-called rational terms resulting from {D — 4) 
times poles in e have to be taken care of for UV singularities, but not for IR singularities 
in this part. Possible rational terms of IR origin are contained in A^si^g which entirely 
consists of 3-point subgraphs and is thus easy to reduce to scalar 2- and 3-point integrals 
-Bo and Cq, thereby keeping the full dependence on D. It turns out that no D-dependent 
prefactors occur in front of IR-singular integrals. In App. |X] we show that this result of 
our specific calculation is not accidental, but generalizes to arbitrary processes at NLO. 
Technically it is easier to evaluate A^sing and A^sing simultaneously according to 
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where A/ = I^^^ — I^'^^ are the differences of the scalar integrals / = Bq, Cq in the two 
IR regularization schemes. Note that IR-finite integrals drop out in AA^^j^ completely. 
Having cancelled UV divergences against counterterms and controlled the D-dimensional 
issues concerning IR singularities, the amplitude is further simplified in four space-time 
dimensions. Specifically, the reduction of SMEs described in App. [B] is performed then. 
The IR-divergent "endpoint part" of the dipole subtraction function, i.e. the contribution 
of the / operator as defined in Ref. [38], is processed through the described algebraic 
manipulations in the same way as LO and one-loop amplitudes. The algebraic Mathe- 
MATICA output of each diagram is automatically processed to Fortran for the numerical 
evaluation. 

Version 2 of the virtual corrections employs FeynArts 3.2 [39] for generating and 
FormCalc 5.2 [40] for preprocessing the amplitudes. The first part of the calculation is 
performed in D dimensions. In particular, the so-called rational terms resulting from the 
UV divergences of tensor loop coefficients are automatically extracted by FormCalc. 
Since the IR divergences that appear in the qq channel are of abelian nature, we exploit 
the fact that they can be regularized as in QED by means of fermion and gauge-boson 
(gluon) masses, rriq and rrig. These masses are treated as infinitesimal quantities (with 
rUg -C rriq) both in the algebraic expressions and in the numerical routines that evaluate 
the tensor integrals, i.e. only the logarithmic dependence on these mass parameters is 
retained. Corresponding IR singularities associated with real emission have been obtained 
from Ref. [38] by means of an appropriate change of regularization scheme. 

Being of diagrammatic nature, the employed techniques are sometimes denoted as 
"brute force" methods. This choice of the terminology might suggest scarce efficiency. In 
fact, the performance of the algorithms is a very important issue that should be assessed 
by means of those quantities that describe the problematic aspects of NLO multi-leg 



calculations: numerical accuracy and CPU time. In this respect our treatment of the 
virtual corrections is characterized by high numerical precision and speed. The numerical 
agreement between the two programs is good, and the CPU time needed to evaluate a 
phase-space point (including sums over colours and polarizations) amounts to about 10^^ 
seconds on a single 3 GHz Intel Xeon processor. This provides a benchmark that can be 
compared with the efficiency of other approaches. 

2.2 Real corrections 

In both evaluations of the real corrections the amplitudes are calculated in the form of 
helicity matrix elements. The singularities for soft or coUinear gluon emission are isolated 
via dipole subtraction [38,41-43] for NLO QCD calculations using the formulation [38] 
for massive quarks. After combining virtual and real corrections, singularities connected 
to collinear configurations in the final state cancel for "collinear-safe" observables auto- 
matically after applying a jet algorithm, singularities connected to collinear initial-state 
splittings are removed via MS QCD factorization by PDF redefinitions. While soft and 
collinear singularities have to be regularized in the "endpoint part" of the subtraction 
function, i.e. the part of the subtraction terms that has to be combined with the virtual 
corrections, no regularization is needed in the subtraction terms for the real corrections. 
In both evaluations the phase-space integration is performed with multichannel Monte 
Carlo generators [44] and adaptive weight optimization similar to the one implemented 
in Racoon WW [45]. 

In version 1 of the real corrections the matrix elements have been calculated using the 
Weyl-van-der-Waerden spinor technique in the formulation of Ref. [46]. Soft and collinear 
singularities are regularized using dimensional regularization. The phase-space integra- 
tion, implemented in C-\ — h, is based on RacoonWW, but the phase-space mappings are 
built up in a more generic way very similar to the approach of LusiFER [47] . 

In version 2 of the real corrections the matrix elements have been generated with 
Madgraph 4.1.33 [48]. As in the corresponding virtual corrections, soft singularities 
are regularized by an infinitesimal gluon mass and collinear singularities by small quark 
masses, which appear only in logarithms in the endpoint part of the subtraction function. 
The Monte Carlo generator is a further development of the one used in COFFER77 [49] 
and for the calculation of the NLO corrections to pp — *> H + 2jets + X [50]. 

In version 2 we have also implemented two-cut-off slicing for the purpose of checking. 
In this approach (as e.g. reviewed in Ref. [51]), phase-space regions where real gluon emis- 
sion contains soft or collinear singularities are defined by the auxiliary cutoff parameters 
^si 5c ^ 1 in the partonic centre-of-mass frame. In real gluon radiation processes, the 
region 
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< fco < 5,^, (2.10) 

where k is the gluon momentum and Vs the partonic centre-of-mass energy, is treated in 
soft approximation. The regions determined by 
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where 9gq is the angle between any hght quark q (including b quarks) and the gluon, 
are evaluated using collinear factorization. We again use an infinitesimal gluon mass and 
quark masses as regulators. In this regularization the contributions of the soft regions for 
light quarks and of the collinear regions can be found in Ref. [52], and the contributions 
of the soft regions involving top quarks can easily be calculated with the explicit results 
for the soft integrals in Refs. [32, 53]. In the remaining phase space no regulators are used. 
When adding all contributions, the dependence on the technical cuts cancels if the cut-off 
parameters are chosen to be small enough so that the soft and collinear approximations 
apply, i.e. the slicing result is correct up to terms of 0{6s) and 0{6c). Since the numerical 
cancellations between the different contributions grow with smaller cut parameters, the 
numerical error blows up if these parameters are too small. 

3 Numerical results 

We consider the process pp -^ ttbb + X at the LHC, i.e. for y/s = 14 TeV. For 
the top-quark mass, renormalized in the on-shell scheme, we take the numerical value 
rrit = 172.6 GeV [54]. All other QCD partons (including b quarks) are treated as mass- 
less particles, and collinear final-state configurations, which give rise to singularities, 
are recombined into IR-safe jets using a fcx-algorithm [55]. Specifically, we adopt the 
fcx-algorithm of Ref. [56] and recombine all final-state b quarks and gluons with pseu- 
dorapidity \ri\ < 5 into jets with separation a/A0^ + Ay'^ > D = 0.8 in the rapidity- 
azimuthal-angle plane. Requiring two b-quark jets, this also avoids collinear singularities 
resulting from the splitting of gluons into (massless) b quarks. Motivated by the search 
for a ttH{H -^ bb) signal at the LHC [24,25], we impose the following additional cuts 
on the transverse momenta, the rapidity, and the invariant mass of the two (recombined) 
b-jets|j pT,h > 20 GeV, \yh\ < 2.5, and mbS > ""^bb.cut- We plot results either as a function 
of mbb,cut oi' foi' "^bb,cut = 0. Note, however, that the jet algorithm and the requirement 
of having two b jets with pT,b > 20 GeV in the final state sets an effective lower limit 
on the invariant mass m^^ of roughly 20 GeV. The outgoing (anti)top quarks are neither 
affected by the jet algorithm nor by phase-space cuts. 

We consistently use the CTEQ6 [57] set of parton distribution functions (PDFs), i.e. 
we take CTEQ6L1 PDFs with a 1-loop running a^ in LO and CTEQ6M PDFs with a 
2-loop running Og in NLO, but the suppressed contribution from b quarks in the initial 
state has been neglected. The number of active flavours is A^f = 5, and the respective 
QCD parameters are Ag^ = 165 MeV and A|^^ = 226 MeV. In the renormalization of the 
strong coupling constant the top-quark loop in the gluon self-energy is subtracted at zero 
momentum. In this scheme the running of Ofg is generated solely by the contributions of 
the light-quark and gluon loops. This yields Q;s(m't)|LO = 0.1178730 . . . and as(mt)|NLO = 
0.1076396 . . .. By default, we set the renormalization and factorization scales, /xr and /ip, 
to the common value fio = nit + m^b cut/2- 



^The experimental analysis of ttH{H —^ bb) will select b quarks with transverse momenta much higher 
than TOb, justifying the approximation mi, = 0. 
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Figure 3: LO cross section (gg + qq) versus contribution from qq annihilation for pp —>■ 
ttbb + X at the LHC as function of the cut 771^,5 cut o'^ the invariant mass of the bb pair. 



3.1 Integrated cross sections 

We first consider results for integrated cross sections. Figure [3] shows the total LO 
cross section and the contribution induced by qq annihilation at the LHC as a function 
of m^b^cut- Fo^^ the chosen setup gg fusion dominates over the qq channel by roughly a 
factor 17. The renormalization and factorization scale dependence of the LO prediction 
is indicated by bands resulting from varying the central scale /iq up and down by a factor 
2 which corresponds to a variation of the cross section by a factor 1.6. Owing to the large 
power of as(/iR)^ in the LO cross section the scale uncertainty is strongly dominated by 
the renormalization scale dependence. We note that the LO cross sections have also been 
reproduced with the program Sherpa [58]. 

Figure H] illustrates the mutual agreement between NLO results obtained with dipole 
subtraction and two-cutoff phase-space slicing. We find that within integration errors the 
slicing results become independent of the cut-offs for 6s ^ 10"^ and Sc ^ 10""^ and agree 
nicely with the result of the subtraction method. While the latter has been obtained 
with 2 X 10* events, the slicing results are based on 10^ events. Still the statistical error 
obtained with the subtraction approach (indicated by the width of the band) is almost 
an order of magnitude smaller than its slicing counterpart (errorbars), demonstrating the 
higher efficiency of dipole subtraction. The results shown in the following are obtained 
with the subtraction approach. 

In Figure Owe show the scale dependence of the LO and NLO cross sections induced 
by the qq channel upon varying the renormalization and factorization scales in a uniform 
or an antipodal way. We observe a sizeable reduction of the scale uncertainty upon going 
from LO to NLO. Varying the scale up and down by a factor 2 changes the cross section 
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Figure 4: Comparison of the relative NLO corrections to pp{qq) — >■ ttbb + X at the LHC 
as obtained with dipole subtraction and two-cutoff phase-space slicing using m^b^cut = 
and /iR = /ip = /io = rrif 
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Figure 5: Dependence of the LO and NLO cross sections of pp{qq) — > ttbb + X at the 
LHC for m^b cut = ^^^ f^o = "^t- 
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Figure 6: LO and NLO cross sections for pp{qq) —>■ ttbb + X at the LHC as function of 
the cut m^b cut o^ the invariant mass of the bb pair, with the bands indicating the scale 
dependence by varying /zr and /ip by factors 1/2 and 2 in a uniform or antipodal way. 

by 55% in LO and by 17% in NLO. At the central scale, the NLO correction is small, 
i.e. ~ 2.5%, and the LO and NLO cross sections are given by ctlo = 85.522(26) fb and 
c"nlo = 87.698(56) fb. The numbers in parentheses are the errors of the Monte Carlo 
integration for 2 x 10*^ events, where the virtual corrections are only calculated for each 
5th event. 

Figure [6] shows the LO and NLO cross sections as function of the cut m^b.cut on the bb 
invariant mass, where the bands indicate the effect from a uniform or antipodal rescaling 
of /iR and /ip by factors 1/2 and 2. The reduction of the scale uncertainty from about 
±50% to ±17% and the smallness of the NLO correction holds true for the considered 
range in wibb,cut; which is motivated by the search for a low-mass Higgs boson. While the 
NLO prediction is consistent with the LO uncertainty band, the shape of the distribution 
is distorted by the corrections. For the central scale we find an NLO correction of ±2.5% 
for small mbb,cut but a correction of —11% for rribbcut = 200 GeV. 

3.2 Differential cross sections 

In this section we consider results for distributions in variables related to the bb pair 
(which in the corresponding signal process pp -^ ttH ± X results from the Higgs decay). 
For each distribution we plot the absolute predictions in LO and in NLO and show the 
relative corrections. These results are based on 2 x 10* events, and no cut on m^^ has 
been applied such that the default scale is Hq = rrit. 

We first show the distribution in the invariant mass m^^ of the bb pair in Figure [3 
The differential cross section drops strongly with increasing m^^, while the relative NLO 
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Figure 7: Distribution in the invariant mass rriy^^ of the bottom-antibottom pair (left) 
and corresponding relative NLO corrections (right) for fi^ = fip = fiQ = m^. 

corrections become large and negative. The increase of the corrections with mbg is larger 
than the one seen in Figure [6] since the scales are fixed to rrtt and not related to m^b- In 
this paper we do not investigate in how far shape distortions induced by the corrections 
could be absorbed into the LO upon using phase-space-dependent scales; we postpone 
this issue until the full NLO corrections including gg fusion are available. The drop of the 
distribution for small ?Tibb is due to the fact that the jet algorithm provides an effective 
cut on this variable. 

The distribution in the transverse momentum Pxbb of ^h^ bottom-antibottom pair 
shown in Figure [S] looks very similar. Again, for our scale choice, the NLO corrections 
reduce the cross section for large values of p-pj^j^. 

Finally, we depict the distribution in the rapidity j/bS of the bottom-antibottom pair 
in Figure [9l In this case, the NLO corrections are rather fiat at the level of 2.5% with a 
slight increase in the backward and forward directions. 

4 Conclusions 

Predictions for the background process pp — > ttbb + X in NLO QCD are indispensable 
for a thorough analysis of ttH production at the LHC 

In this paper we have made the first step towards the full NLO calculation upon evalu- 
ating the contribution from quark-antiquark annihilation. We made use of the Feynman- 
diagrammatic approach augmented by recently developed reduction techniques for one- 
loop tensor integrals. We have devoted an appendix to the general issue of rational terms 
resulting from ultraviolet or infrared (soft or collinear) singularities within dimensional 
regularization. In particular, we have shown that rational terms of infrared origin cancel in 



11 




^-1[%] ppiqq) ^ tthh + X 




50 100 150 200 250 300 350 400 
PT,bb [GeV] 



50 100 150 200 250 300 350 400 

PT,bb [GeV] 



Figure 8: Distribution in the transverse momentum p^ ^^^ of the bottom-antibottom pair 
(left) and corresponding relative corrections (right) for /xr = /xp = /^o = ""^t- 
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Figure 9: Distribution in the rapidity y^^ of the bottom-antibottom pair (left) and cor- 
responding relative corrections (right) for /ir = /xp = /^o = "^t- 
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truncated one-loop diagrams for arbitrary processes in the Feynman gauge and thus result 
only from wave-function renormalization. Based on this observation we have formulated 
a general recipe for the determination of rational terms in one-loop amplitudes. 

Our calculation demonstrates that the Feynman-diagrammatic approach can be suc- 
cessfully applied in the context of six-particle processes at the LHC, providing excellent 
numerical stability and high speed. The CPU time needed to evaluate the full virtual 
corrections with a single processor is of the order of 10~^ seconds per phase-space point. 
Based on these encouraging results, we expect to be able to extend this calculation to the 
technically more challenging gluon-fusion channel. 
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Appendix 

A Rational terms in one-loop amplitudes 

In this appendix we elaborate on the issue of so-called rational terms that occur 
within dimensional regularization when algebraic factors depending on the space-time 
dimensionality D multiply loop integrals that contain UV or IR (soft and/or collinear) 
singularities, which give rise to poles in D—A. The calculation of rational terms is of partic- 
ular importance for so-called unitarity or generalized-unitarity methods (see Refs. [10-13] 
and references therein). In this context, these terms are typically obtained by recursion 
relations [59-61] or by exploiting the full D-dimensional dependence of the tree ampli- 
tudes [12,62,63]. Explicit recipes to derive rational terms in the context of specific 
methods, which make use of loop integrals in shifted space-time dimensions or employ a 
numerical reduction at the integrand level, can be found in Refs. [64,65]. Here we discuss 
rational terms in the framework of one-loop calculations that employ an arbitrary set of 
tensor (and scalar) loop integrals in D = 4 — 2e dimensions and derive general properties 
that are independent of the explicit algorithm used for tensor reduction. 

Specifically we investigate rational terms resulting from unrenormalized truncated loop 
amplitudes and do not consider external self-energy corrections (wave-function renormal- 
ization constants). Since the latter enter via derivatives, our arguments cannot be applied, 
but these contributions are easily calculated once and for all. We classify the different sit- 
uations in which rational terms arise and describe simple procedures and results for their 
actual calculation. In particular, we demonstrate that — for any scattering amplitude in- 
volving quarks and gluons — the rational terms originating from IR poles cancel within 
individual Feynman diagrams. This important property implies that, after separating the 
rational terms of UV type, the coefficients of all IR-divergent tensor A^-point integrals 
can be evaluated in four dimensions. In practice the wave-function renormalization con- 
stants represent the only source of rational terms of IR origin. This greatly simplifies the 
algebraic manipulation of IR-divergent scattering amplitudes. 
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A.l Classification of rational terms 

Algebraic factors containing the dimensionality D result from two different sources in 
one-loop amplitudes: 

1. "Trace-like" contractions among metric tensors or with Dirac structures lead to 
expressions such as g^^g^^ = D, 'j'^^l'^d^Li' = l^'I'-ln = (2 — D)(/i, etc. 

2. In the reduction of tensor one- loop integrals to standard scalar integrals (such as 
the usual Passarino-Veltman reduction [35]) the tensor coefficients are eventually 
obtained as linear combinations of the scalar integrals which form a basis of func- 
tions. In such linear combinations, the tensor coefficients containing metric tensors 
in their corresponding covariants receive prefactors with a dependence on D. 

Thus, we can distinguish four different types of rational terms, classified according to 
type 1 or 2 being of UV or IR origin. 

We employ the notation of Ref. [5], where the covariant coefficients of A-point integrals 
with rank R are denoted as Tj^ ,• . It is convenient to treat the UV- and IR-divergent 
parts of tensor integrals separately. To this end, we write 

Tl" i =t^ i ^flii^AiL (A.l) 

euv 

where -R^ j represents the (IR-finite) residue of the UV pole of T^^ j , and T^^ ^ is free 
from UV divergences but can contain single and double poles in ejR resulting from soft or 
collinear divergences. The only IR-divergent 2-point functions are those without a scale. 
These represent a special case since they vanish as a result of cancellations between IR 
and UV poles, i.e. they are formally proportional to (1/euv — I/cir). In order to separate 
these UV and IR poles, we isolate the UV divergences by writing, in the notation of 
Ref. [5], 

5o...o 1...1 (0, 0, 0) = 5o...o 1...1 (0, 0, 0) + "^-l^^r^. (A.2) 



n + 1 euv 



The UV-subtracted part is IR divergent. 



-lY 5,. 



mO 



5o...oi...i (0,0,0) = -^^^, (A.3) 



n 



1 eiR ' 



and exactly cancels against the UV pole. However, we do no set scaleless 2-point integrals 
to zero and treat the rational terms resulting from 1/euv and -B...(0, 0, 0) separately. Since 
scaleless 2-point functions require light-like momentum transfer (p^ = 0), such integrals 
only occur in external self-energy corrections, i.e. in wave-function renormalization con- 
stants, and in the reduction of higher-point functions {N > 2). In the latter case, as we 
will show below, the rational terms of IR origin cancel out. 

A.2 Rational terms of UV origin 

The residues Rf-^,,,i^ of the UV poles of general one- loop tensor integrals are simple 
polynomials of the external momenta and their explicit form is well known (see, e.g., 
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Figure 10: Kinematical configurations for soft (left) and coUinear (riglit) IR singularities 
in one-loop diagrams. 

App. C of Ref. [4] and App. A of Ref. [5]). In particular, in contrast to IR divergences, 
the UV poles do not depend on kinematical properties of the amplitude such as on-shell 
relations of momenta. This renders rational terms of UV origin very simple: Terms of 
type 2 can be directly included in the tensor reduction in a generic way, as e.g. done in 
Refs. [4,5,35,53], and terms of type 1 can be extracted during the algebraic reduction of 
each Feynman diagram by means of a trivial expansion. 



N 



fiD)Tl,,^ = f{D)Tl,,^ + 



m 

euv 



- 2/'(4) 



R 



N 



(A.4) 



In the following we discuss the remaining rational terms that result from f{D)T^ ^ when 



T!^ ,„ contains poles of IR origin. 



-11. ..In 



A. 3 Rational terms of IR origin 

IR divergences of one-loop integrals are more complicated than UV singularities, since 
they depend on specific kinematical properties of amplitudes. According to Kinoshita [66] 
they can be classified into soft and collinear singularities as indicated in Figure [TOl 

• A soft singularity arises if a massless particle is exchanged between two on-shell 
particles (see l.h.s. of Figure [TOl) . The singularity is logarithmic and originates 
from the region in momentum space where the momentum transfer of the massless 
propagator tends to zero. 

• A collinear singularity arises if an external line with a light-like momentum (e.g. a 
massless external on-shell particle) is attached to two massless propagators (see r.h.s. 
of Figure [TOl) . The singularity is also logarithmic and originates from the region 
in momentum space where the loop momentum of the two massless propagators 
becomes collinear to the momentum p of the external particle. 

We first consider rational terms of IR origin that can result from the reduction of 
tensor integrals (type 2). As can, for instance, be seen from the results of Ref. [5], only 
tensor coefficients T^ whose covariants involve the metric tensor can get D-dependent 
prefactors in reduction identities and could, thus, lead to rational terms of type 2. In 
Section 5.8 of Ref. [5] it was, however, shown that these tensor coefficients are IR finite. 
Thus, no type 2 rational terms of IR origin can result at all. 
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Moreover, reparametrizations of tensor integrals resulting from shifts of the loop mo- 
mentum and permutations of the propagators do not give rise to /^-dependent factors 
or relations between (IR-finite) tensor coefficients of type T^ and IR-singular integrals. 
Therefore, in order to demonstrate that a certain diagram is free from IR rational terms 
it is sufficient to find for each soft- or coUinear-singular region a specific representation 
that is manifestly free from IR rational terms of type 1, i.e. to find an expression in which 
the corresponding IR-divergent part is expressed as a linear combination of IR-divergent 
tensor (or scalar) integrals with coefficients that are independent of Z^. In explicit cal- 
culations, this can be achieved by means of an algebraic reduction that implements all 
possible relations between the IR-divergent parts of tensor integrals. This task is non- 
trivial since the standard scalar integrals do not provide a unique representation of IR 
divergences. For instance, IR-finite parts of a diagram can be expressed as linear combi- 
nations of IR-divergent 4-point and 3-point scalar integrals. More generally, IR-singular 
A^-point tensor integrals can be expressed in terms of IR-divergent 3- and 2-point scalar 
integrals plus IR-finite terms. This reduction of IR singularities is explicitly implemented 
in the algorithm presented in Ref. [37] and can be summarized by the following formula 
(see Eq. (3.14) in Ref. [37]) which relates the IR-divergent part of iV-point tensor integrals 
to 3-point tensor integrals associated with the IR-divergent triangle subdiagrams, 

Tl^i...t,RiPo, • • • , PN-umo, . . . , tun-i) 

N-l N-1 

= Yj Yl AnkC^^,„^^{pn,Pn+i,Pk,mn,mn+i,mk) + IR-&mtepaTt. (A. 5) 

n=0 '==0 

fc 7^ 71 , ri + 1 

The coefficients Ank are independent of D since all relations between IR-divergent tensor 
integrals are free from D-dimensional coefficients. For the explicit form of Ank and details 
of the notation we refer to Ref. [37]. In practice, using (lA.Sp and performing a subsequent 
reduction to scalar integrals, one can construct a unique representation of IR divergences 
in terms of 3- and 2-point scalar functions. The 3-point functions on the right-hand 
side of flA.Sp can be subtracted from the complete tensor integral leading to an IR-finite 
expression which can be evaluated in 4 dimensions. Thus, only the re-added 3-point 
functions flA.Sp have to be manipulated in D dimensions. 

Using this approach, we have observed in explicit calculations for pp -^ ttH [27], 
pp — i> tt + jet [28], and pp — > WW + jet [16] that — after complete reduction of the IR 
divergences — the coefficients of the IR-singular Cq and Bq functions are independent of 
D, i.e. that rational terms of IR type cancel completely. This indicates, a posteriori, that 
the terms f{D)T^^^^ in (lA.4p can be replaced by f{4:)Tj^^^-^ from the beginning in the 
calculation. In the following we demonstrate that, in the Feynman gauge, the cancellation 
of IR rational terms is a general property of scattering amplitudes involving an arbitrary 
number of external quarks and gluons. 

To this end, we inspect the integrand of a general one-loop IR-divergent diagram in 
momentum space and, in the spirit of Ref. [37], we separate the IR singularities associated 
with different soft and collinear regions and relate them to triangle subdiagrams. Using 
on-shell relations we show that, in the soft and collinear regions, the integrands can be 
cast into a form that is free from "trace-like" contractions, which potentially produce In- 
dependent factors. In this way we obtain a generic representation of the IR singularities 
that is manifestly free from rational terms of IR type. 
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Figure 11: Subdiagrams containing soft-singular integrals, with solid lines indicating 
quarks, epicycles gluons, and dotted lines ghosts. 



The considerations presented in the following are not new. A similar approach is 
used, for instance, in Ref. [67], where IR and UV singularities are subtracted at the 
diagrammatic level and isolated in simple process-independent terms in order to obtain 
numerically integrable expressions. Using similar techniques, here we consider soft and 
coUinear contributions from the perspective of an analytic calculation in terms of divergent 
one-loop tensor and scalar integrals and discuss the rational terms associated with IR 
divergences. 

(i) Soft singularities: 

Figure [TT] shows all potentially soft-divergent subdiagrams containing quarks, gluons, 
or Faddeev-Popov ghosts. Since the soft singularity is related to zero-momentum transfer 
on the internal line linking the two external particles, it is convenient to identify the 
integration momentum q of the loop integral with the momentum on this line. Since the 
soft singularity is logarithmic, all contributions of q in the numerator of the integral are 
IR finite. In other words, being only after the IR-divergent part we can set q to zero in the 
numerator and in propagators that do not belong to the soft-singular triangle subdiagram. 
This procedure immediately kills the three subintegrals (b)-(d) with a quark on the q line 
because of the factor ^ in the (massless) quark propagator. Diagram (g) with a ghost 
coupling to external gluons does not contribute either, because the integrand receives 
factors (depending on the gauge) of qSa — * or paEa = from the coupling of the ghost to 
the on-shell gluon with momentum pa and polarization vector £„. The amplitude A1(8a) of 
diagram (a), of course, involves a soft singularity, but without any D dependence, as can 
be seen in the following example where we consider an outgoing quark-antiquark pair. 



M 



(Sa) 



+ :^a + ma 



/Q 
d^q^^Ua{Pa)lt, , , N2 2 

q^ [q + Pa) -ml 

d q^- Ua{Pa) 7m 
J u 



{q + Pa) 



r(g) 
r(o) 



{q-pbY -ml 

-^b + mb 
{(l-PbY -ml 



luVb{pb) 
luVbiPb) + ••• 
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The Dirac structure r(g) contains the remaining part of the diagram and the elhpses 
stand for terms that are not singular if the gluon becomes soft. For the last equality the 
Dirac equation was used twice. The soft singularity, which is just contained in the scalar 
3-point function, does not receive D-dependent factors and, thus, does not deliver rational 
terms, because r(0) is a tree-like structure and does not contain a trace-like contraction 
that would lead to factors of D. Such factors, e.g., arise in terms like 7'^r(0)7^, which 
are absent in the soft-singular part. The same reasoning applies also to all other possible 
fermion-number flows in diagram (a). The remaining two diagrams (e) and (f) of Figure [10] 
can be analysed in the same way, leading to the same conclusion that no D-dependent 
factors multiply soft-singular integrals. For brevity we show this only for diagram (f) 
explicitly, 

.. f.D 9^" ^f,* a^ujq + 2pa)A - guxi'^q + pg)^ + gx^jq - Pa)u -p^^^ ^ 
M(8i) = / d q ^^e'^ — ■ -^ 1 {q) 

' ' J q^ [q + PaY 

^ p*9pr{-q+^^Pb)a + 9ra{2q-pb)p-gap{q+Pb)r 



{Q-PbY 



J-D 9 n^'^gpi^PajX 9y\Pa,p 9X/^Pa,u -pi^r fri\ 

q^ [q + PaY 

p* '2gprPb,a - QraPb^p " 9cTpPb,- 

{q-PbY 






.D ['^iPaPb)eluelr - '^{CPh)Pa,uelr " HPae*b)el^Pb,r + {elel)pa,^pb,r]T''^{0) 
d q ■ '■ ■ ■ 

q'^{q + PaYiq-PbY 

+..., (A.7) 

where we assume the external on-shell gluons to be outgoing. The last line of this result 
cannot contain explicit factors of D, since those would require a trace-like contraction 
r'^,^(0); instead F is only contracted with momenta pa,b and polarization vectors e*j,. 
Formally the contraction r'^y(O) occurs, but with a proportionahty to {pael){pbel) which 
vanishes owing to the on-shell condition of the gluons. 

It is well known that the soft singularities are ruled by the eikonal current, with the re- 
sult that divergences connected to soft-particle exchange between a and b are proportional 
to (paPb)- For ^A(sa) this factor is already explicit in flA.6p . for A^(8f) in flA.7l) this fact 



becomes obvious after making use of the Ward identities pa^u^l r^'^'^ i^) = ^l uPb,T^'^^ {^) = 
Pa,uPb,T^'^^ (0) = 0, which are valid if all other external particles are on shell and all 
diagrams contributing to F are summed over. 

(a) Collinear singularities: 

Figure [12] shows all potentially coUinear-divergent subdiagrams containing (massless) 
quarks, gluons, or Faddeev-Popov ghosts. We identify the integration momentum q with 
the momentum of one of the two propagators attached to the external light-like line, which 
carries the momentum p {p^ = 0). The collinear singularity stems from the region in q 
space where q is collinear to p, i.e. q^ = x{q)p^ with some scalar function x{q). Since the 
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Figure 12: Subdiagrams containing collinear-singular integrals, with solid lines indicating 
quarks, epicycles gluons, and dotted lines ghosts. 

singularity is only logarithmic, we do not change the singularity if we replace g^ by x{q)p^ 
in the numerator of the amplitude. The on-shell conditions (Dirac equation for quarks 
and transversality condition for gluons) then imply that the collinear-singular parts of 
diagrams (a)-(d) in Figure [12] do not receive D-dependent factors: 



jj 2[1 + x{q)\u{p)p^T^{x{q)p) 



,. - ^v-^7Trt5- '"<"<«>"' + ■■■ 



+ ... 



M 



(9b) 



J q^iq + pY 



(A.8) 



M 



(9c) 



+ ... 



(A.9) 



D„ ^M* 3f^v{(l + 2p)a - 9v\{'^q + P)^, + 9xM - p )^ ^^a 

q^{q + pY 



T''\q) 



= I d^g e^- 

^Ja^q 

= I d^g e^- 
= + 



. 9f, u{x{q) + 2)pa - ^^A(2x(g) + l)p^ + gxf,{x{q) - l)p^ .p^;,^ ( ^ ^ , 

Q' (Q' + P) 
[(a;(g) + 2)£:pa + (x(g) - l)^!^.] r-^(x(g)p) 



9m 



g2(g + p)2 
. a;(g)p^ 
g2(g + p)2 



g2(g + p)2 

r(g) 

r(x(g)p) 



+ 



(A.IO) 



(A.ll) 



Note that diagrams (b) and (d) do not have collinear singularities at all. The collinear 
divergences in diagrams (a) and (c) do not receive Z)-dependent factors since, in the 
collinear region, the F terms do not take part in trace-like contractions like 7^F'^ or 
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V^u, and terms q^^ . . . q^" inside F yield tensor structures [x{qy\"'p^^ . . .p^" without metric 
tensors. The above IR-divergent integrals are easily expressed in terms of usual 3-point 
functions by observing that in the collinear limit all propagators inside T{x{q)p) are linear 
in x{q) owing to {q + PkY — ^ 2a;(g)ppfc + 'Pk- Thus, one can easily express the denominator 
of T{x{q)p) as a linear combination of single propagators via partial fractioning, as done 
in Ref. [37]. 

In summary we have obtained a representation of the soft and collinear singularities of 
generic Feynman diagrams in terms of 3-point tensor integrals that are explicitly free from 
D-dimensional prefactors. We can, thus, conclude that rational terms of IR origin cancel 
in any unrenormalized scattering amplitude and can be neglected a priori in explicit cal- 
culations. This property is a consequence of the Lorentz structure of the gluon couplings 
and the logarithmic nature of IR singularities within the conventional Feynman gauge. 
However, the cancellation of IR rational terms holds also in more general gauge-fixings, 
as for instance the background-field Feynman gauge (see Ref. [68] and references therein), 
where the Lorentz structures in the gluon couplings and the poles of the propagators 
behave as in the Feynman gauge. 

We point out that the cancellation of IR rational terms is independent of the actual 
reduction method employed in the calculation and is generally valid in any approach 
where the IR-divergent parts of loop diagrams are entirely expressed in terms of tensor 
(or scalar) A^-point integrals in 4 — 2e dimensions. 

A. 4 A recipe for determining rational terms 

Based on the above considerations we can formulate a simple algorithm for determining 
all rational terms of either UV or IR origin. For unrenormalized truncated one-loop 
diagrams, i.e. excluding counterterm diagrams and self-energy corrections to external 
lines, proceed as follows: 

1. Separate UV and IR divergences of all tensor integrals as in flA.ip . thereby keeping 
track of the poles of all scaleless 2-point integrals as indicated in flA.2p and flA.3p . 



2. Extract the rational terms of UV origin as described in flA.4p including the rational 



terms resulting from UV poles of scaleless 2-point integrals. 

3. Ignore rational terms of IR origin upon replacing f{D) -^ /(4) on the right-hand 
side of (lA.4p . because the arguments of the previous section show that all rational 
terms resulting from [f{D) — /(4)]T^^ ■ in (lA.4p compensate each other (even if 
they may arise in intermediate steps). 

This recipe does not apply to wave-function renormalization constants. These can be 
easily calculated, and exphcit results are, e.g., given in Eqs. (2.27)-(2.28) of Ref. [27]. 

B Four-dimensional reduction of Dirac chains to standard matrix elements 

Here we outline the algebraic procedure employed to reduce Dirac structures to stan- 
dard matrix elements. The reduction is based on the strategy described in Sect. 3.3 
of Ref. [29], which is worked out for massless six-fermion processes, and involves addi- 
tional features to treat the Dirac chains associated with massive top quarks. This method 
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exploits a set of relations that do not give rise to denominators involving kinematical vari- 
ables, thereby avoiding possible numerical instabilities in exceptional phase-space regions. 
Refraining from a detailed description of the entire reduction algorithm, which is quite 
involved, we only outline the basic principles, which can be traced back to a few simple 
identities. 

Tree and loop diagrams give rise to a large number of Dirac structures of the type 

vipih'"' • • .#n • ■■u{p2) v{p3)Y' ■ ■■i>h ■ --uiPi) v{p5)Y' • • -^fci • --uipe), (B.l) 

which consist of gamma matrices (or slashed momenta) that are sandwiched between 
the spinors v{pa) and u{pa+i) of the six (anti)fermions. For convenience we consider 
the crossed process ggttbb -^ 0, where all particles and their momenta are incoming. 
While the chains associated with massless quarks (a = 1, 5) contain only odd numbers of 
Dirac matrices, inside the top chain (a = 3) also even numbers of Dirac matrices appear. 
The open Lorentz indices /ij , Uj , pk in (IB. II) are always pairwise contracted via the metric 
tensor. 

The combinations of Dirac matrices occurring inside individual chains v{pa) ■ ■ ■ u{pa+i) 
can easily be simplified by means of elementary relations: 

(i) Standard Dirac algebra permits to reduce 7^ ... 7^ contractions inside Dirac chains, 
to bring gamma matrices and ^ terms into a standard order through anti-commu- 
tations, and to eliminate ^^ terms. 

(ii) All ^a and ^a+i terms can be eliminated with the Dirac equation. 

(iii) The ^i terms associated with one of the other external momenta {i ^ a,a + 1) can 
be eliminated via momentum conservation. 

After these simplifications, which we perform in D dimensions, we are left with a 
still large number of Dirac structures of 0(10^). To obtain a further reduction we employ 
additional identities which permit to shift ^ terms and 'j^ matrices with open indices from 
one chain to another, thereby permitting further simplifications of type (i)-(iii). This part 
of the reduction relies on four- dimensional relations and is performed after separation of 
all {D — 4)-poles in dimensional regularization. 

All four-dimensional identities are derived from basic relations which follow from 
Chisholm's identity (see Eqs. (3.4)-(3.6) in Ref. [29]) and read 

7'^7"7'^u;± ® 7^ = -f^'uj^ ® (7^7^7"^± + 7"7^7^a;^) , 

7"7'^7'^cu± ® 7^ = Y^± ® (7'^7/.7"^± + 7"7/.7^^t) > 

7"7^7^a;± ® 7m = l''^± ® (7''7°7/.^± + 7m7°7^^t) , (B.2) 

where u± = (1 ± 7^)/2 are the chirality projectors and the tensor products connect 
different Dirac chains. In order to exploit these relations we introduce chirality projectors 
inside every Dirac chain, v{pa)Tu{pa+i) = J2\=±^{Pa)^^\u{pa+i)- Then we can use flB.2p 
to exchange 7^7^-terms between chains that are connected via 7^-contractions. A simple 
application of (1B.2P is given by 

7'^7"7''u;± ® 7^^7^7m = 4^°^7''^± ® 7m^t + ^7'^^± ® 7"^±- (B.3) 
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These identities permit to eliminate double Lorentz contractions between two Dirac chains. 
Alternatively we can use (1B.3P in the opposite direction, in combination with the Dirac 
equation. This yields the relations 

- -Yl^'ib^T^iPa) ® -f^,uj^u{ph) - —1^uj±u{pa) ® lja^±u{pb) + m^mb-tcrm, 

V{pa)'^b^± ® v(pb):^a'^T = {P aPb)v {p a)l^ ^ ± ® v{pb)-f^UJ^ 

+ -Y'^{pa)ibl^^± ® v{pb)-f^,uj^ + —v{pa)'j''uj± ® v{pb)'^ al,,^^ + m^mb-term, 

V{pa)'^b^±®i>a^±u{pb) = {PaPb)v{Pa)l^^±®l,j,^±u{pb) 
Tfl Tflh 

+ -Yv{pa)ibl^UJ± ® Jf,UJ±u{pb) + —v{j)a)l^UJ± ® :^alfiUJ^u{pb) , (B.4) 

where the terms proportional to marrib vanish since the (anti)spinors associated with 
particles a and b belong to different Dirac chains and, thus, at least one of them is 
massless in our case. The relations (IB. 40 can be used to reduce the number of '^i terms 
in the Dirac chains. 

We give two explicit examples to illustrate the four- dimensional reduction of Dirac 
structures that involve one massive Dirac chain: 

Example 1: 

We reduce terms involving double contractions of the type ''y^Y®ln®lv to structures 
involving only single contractions of Lorentz indices as much as possible. This 
procedure somewhat generalizes Step 1 in Section 3.3 of Ref. [29]. Following the 
notation of that reference, we use the shorthand [F]^. = v{pi)Vujpu{pj) and consider 
Dirac structures of the type 

where the terms Ai consist of ^^-products, each containing rti slashed momenta. By 
means of the following two steps the structures flB.5l) can be recursively reduced 
to [7^7^^] 34 [7^] 12 [7^^156 ^^'^ terms that are free from double contractions such as 

Step 1: \i rii > 1 ioi i = 1 or 2, then we write Ai = Aij^alfb and, using (IB. 20 . we 
shift '^ai'b to the chain that contains Aq'-j'^'j". Then we perform the simplifications 
(i)-(iii) in four dimensions. This step is iterated until ni,n2 < 1. 

Step 2: If ni, ^2 < 1 and uq > 0, then we can write Aq = Ao^a with a G {k, I, m, n}, 
since ^jj are eliminated by means of (i)-(iii), also in four dimensions. In this case, 
using (IB. 20 . we shift ^a and one of the matrices 7^,7"^ from the Ao-chain to that 
y4j-chain where we can eliminate "^a by means of the Dirac equation and other 
simplifications (i)-(iii). Then we restart with step 1. 

This procedure recursively reduces the number of |^-terms no+^1+^2 until rii, n2 < 1 
and no = 0, which automatically implies uq = rii = n2 = since only one of the 
three Dirac chains (the massive one) can contain an even number of Dirac matrices. 
Thus the only 7^7"^ ® 7/^ ® 7j/ structure that survives is [7^7^^] 34 [7/i]i2 b^lse- 
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Example 2: 

We consider Dirac structures of the type [7'^]^2 [7M^fc]34 [Alse- Using the relations 

b']% bM. = +^ blf^ ^21^41^4 = +^ b'hMu b£. = o, (b.6) 



which follow from ( ]B.2p . and using momentum conservation, we can achieve that 
the index k takes only the values k = 5,6 for each chirality configuration {par). 
Eliminating one pi via momentum conservation, the index / can take three values, 
leading to six different index pairs {kl) per chirality configuration. Two out of the 
six possibilities can be easily eliminated by relations like ( ]B.4p : 



Ms] 34 lAhe = (P4P5) [7^7^.134 bise - -Y bf^li^hhi b'^' 

b^^Mu IMle = (PiPe) [7m7^]^4 blJe " ^ b^^'y^^Mu b'^ 



b^^h]M [h]% = -iP3P5) [7^7/.] ^4 [71^6 - -J- [|^57^7m]^4 bite + 2^5,^ [1]^4 [h]% ' 

[7mA]^4 [h]l6 = -(PsPe) [7^7/.]^4 [71^6 " ^ [^7^7m]^4 bl^e + 2^6,^ [1]^4 [h]56 ' 

(B.7) 



where the use of (1B.4P in the last two equations required an anticommutation of 



7/i 1^5,6 leading to additional contributions, and the terms proportional to m3_4 on 
the right-hand side of (]B.7|) can be further reduced as in Example 1. Another (kl) 



combination can be eliminated by using identities like (lB.4p for the chains [. . .]i2 
and [. . .jse after shifting i^^ to [. . .]i2. In order to achieve this, one has to apply the 
Dirac equation for the massive fermion inversely in the first step: 

b,]t2 b'Mu [^1]% = ^ b,]% b'hMl 1^]% = ^ IMa.t b']l [^i]% 

= ^ [7.^7.]f, [7IJ4 [7l£ = ^ b.]% b'lMl b']% 

= (PiPe) [7M]f2 b^lujti bl% ' 
b,]t2 b''h]u Me = ■■■ = (P1P5) b,]t2 [7^.]^4 bll , 

b£2 b'Ml Me = ■■■ = iP2P6) b£2 [7^.1^4 bll , 

b.t b''h]l Me = ■■■ = (P2P5) b.t [7^.] J4 [71^6 • (B.8) 

One additional relation per chirality configuration results upon exploiting = 
[7'^(^ + ... + MY]i2bf^Mub'^]le similar to Step 5 in Section 3.3 of Ref. [29], 
however, this procedure is quite tedious. 

The complete reduction algorithm consists of several procedures of this type, each con- 
sisting of combinations of the identities (]B.2p - (lB.4l) and the operations (i)-(iii). 

In the case of massless 6- fermion processes [29], all Dirac structures were reduced to 
10 types of SMEs of the form 
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Counting the different chiralities p,a,T = ±, which yield 8 or less combinations per type 
of SME depending on the type, the total number of independent "massless" SMEs was 
80. In addition to these SMEs, the qq -^ ttbb reduction yields 15 types of SMEs of the 
form 

b%bMMu^ b%b%ib,i.]u, ^aY^.b^bMu, W&li^^lIJilL, (b-io) 

where the chain [. . .]34, i.e. the top-quark chain, involves an even number (0 or 2) of 
Dirac matrices. In the two independent reduction algorithms that we have implemented 
the total number of SMEs for qq -^ ttbb, counting all types flB.9l) - flB.10p and chiralities 
p, 0", r, is 148 and 156. As it is obvious, the presence of the top mass increases the number 
of the SMEs by roughly a factor 2. Moreover, also the complexity of the form factors 
associated with each SME grows considerably with respect to the case where all fermions 
are massless. 

C Benchmark numbers for the virtual corrections 

In order to facilitate a comparison to our calculation, in this appendix we provide 
explicit numbers on the squared LO amplitude and the corresponding virtual correction 
for a single non-exceptional phase-space point. The set of momenta for the partonic 
reaction qq — > ttbb is chosen as 

Pg = (500,0,0,500), 

p'^ = (500,0,0,-500), 

p(^ = (327.5045589027869,107.1276753641986,-107.9290580423663,-233.1168284428635), 

Pi = (276.6425142763093,-107.4949148022111,153.8289259355409,-107.3397668261919), 

pi^ = (233.9459027189062,82.55875671042013,-77.70592645955253,204.6375480757531), 

Pg = (161.9070241019976,-82.19151727240762,31.80605856637796,135.8190471933023), (C.l) 

with the components given in GeV and rrit = 172.6 GeV. We give numbers on the spin- 
and colour- averaged squared LO amplitude jA^'^^p and on the sum of the relative virtual 
NLO correction ^virt and the contribution 6i of the / operator of the dipole subtraction 
function as defined in Ref. [38]. In more detail, we split the relative correction into a 
contribution originating from closed fermion loops, 5ferm (comprising contributions from 
the gluon self-energy, the triple-gluon vertex correction, and the renormalization constant 
of the strong coupling), and the remaining loop corrections, called ^bos, and 6i. Note that 
for qq — > ttbb the fermionic part is IR finite, while ^bos is IR divergent. Adding 6i to 
(5bos or (5virt = (^ferm + ^bos, all IR divergences cancel, and the sum is independent of the 
IR regularization scheme. The values of the strong coupling constant at /xr = rrit in the 
setup described in Section [3] are 

tts(mt)|Lo = 0.1178730139006150, as(mt)|NLO = 0.1076396017050965. (C.2) 

At the phase-space point (JC.ip we find 

|>lLO|2/^8 _ 0.4487410759198035 • 10-^ GeV-^ 
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Ovirt+I 
Ovirt+I 
Oferin 
"ferin 
^bos+I 
^bos+I 



= 0.4487410759198011 ■ 10"^ GeV"^ 

Madgraph 



= -0.1290522911043483, 

vcrsionl 

= -0.1290522911137204, 

version2 

= -0.06326213639716407, 

versionl 

= -0.06326213639715421, 

version2 

= -0.06579015470718421, 

vcrsionl 

= -0.06579015471656619, (C.3) 

version2 



where we divided out the strong couphng constant Qs from jA^^'-'p. The agreement 
between our two independent versions of the virtual corrections is typically about 10 digits 
at regular phase-space points. 
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